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Abstract 

Solution of the Cox-Thompson inverse scattering problem at fixed energy [T] [2l [3] is reformulated 
resulting in semi-analytic equations. The new set of equations for the normalization constants and the 
nonphysical (shifted) angular momenta are free of matrix inversion operations. This simplification is a 
result of treating only the input phase shifts of partial waves of a given parity. Therefore, the proposed 
, method can be applied for identical particle scattering of the bosonic type (or for certain cases of 

identical fermionic scattering). The new formulae are expected to be numerically more efficient than 
the previous ones. Based on the semi-analytic equations an approximate method is proposed for the 
generic inverse scattering problem, when partial waves of arbitrary parity are considered. 
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1 Introduction 



Inverse quantum scattering theories represent an important field of physics research for more than fifty 
' years. The reason for this ever growing interest is simple: application of the inverse methods may yield 

information about the interactions governing the microscopic processes which can be detected by studying 
simple macroscopic scattering experiments. 

There are several theories of inverse scattering; in this paper we treat only one of them, namely the 
Cox-Thompson (CT) method at fixed energy[T]. This procedure produces potentials that possess finite 
values at the origin and non-zero first momenta. Furthermore the CT scheme for inverse scattering at 
fixed energy proves to be a very efficient constructive procedure ever proposed, although it becomes 
difficult to carry it through if input data possess large errors. Thus simplifications of this method is a 
question of interest. 

In this paper we show how the CT method can be simplified when it is applied to identical particle 
scattering of bosonic type. In this case only the even partial wave {£ — 0,2,...) phase shifts can be 
extracted from the experimental data and the solution of the nonlinear CT equations can be simplified. 
The simplification means here that while the general nonlinear system of the CT equations contains 
implicit matrix inversions, the simplified equations do not. In the case of even partial waves (and also for 
the odd waves, ^ = 1, 3, ...) the matrix inversions can be performed analytically. Therefore the simplified 
method is called as 'semi-analytical' one. 



1 



Explicit calculations show that the new method provides more stable solutions to the nonlinear CT 
equations. Also, it turns out that the separate treatment of odd and even partial waves and the subsequent 
addition of the resulted two potentials is a good approximation. 

In the next section we review the CT method. The simplified equations valid in case of even or odd 
angular momenta are derived in section 3 and various applications are presented in section 4. Section 5 
is left for a short summary. 



2 Cox-Thompson method 

We consider the scattering of particles by a spherically symmetric potential. The main equation of our 
interest is thus the radial Schrodinger equation. The well-known partial wave decomposition leads to 

for the ipi{r) partial waves involved in the scattering event. 

Let us introduce the commonly used dimensionless quantities, x = kr for the distance, and q(x) = 
^^'^/'^^ for the potential where E = ffk'^/2m denotes the scattering energy, m the reduced mass, and k 
the wave number. 

One can rewrite the radial Schrodinger equation into a more appealing form 

Dg{x)M^) ^Ki + '^)Mx), les, (2) 

where the differential operator Dq{x) is defined as 



..^ (. 
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D,ix) ^ x' + 1 - qix) j . (3) 

Clearly in case of a zero potential the two linearly independent solutions of the Schrodinger equation ([2]) 
are the Riccati-Bessel and the Weber-Schlafli functions defined as 

= y^J/+i(x), vi{x) ^ ^J^Yl^l{x). (4) 

It is a well-known fact that one can transform the free solutions ui (x) into the solutions V-"; {^) by using 
a generalized translation operation [S]: 

ij;i{x)=ui{x)- I dtt-^K{x,t)ui{t), I € S. (5) 



Jo 

This equation is often referred to as the Povzner-Levitan representation of the wave function, the quantity 
K{x,y) is named as the transformation kernel. 

The inverse potential is calculated from the transformation kernel as follows 

2d K{x,x) 

The transformation kernel is determined by an integral equation of the Gcl'fand-Levitan-type (GL-type) 
which contains the symmetric input kernel 3(a;,y]j: 



K(x,y) ^ g{x,v) - dtr^ K(x,t)g{t,y), x>y. (7) 
Jo 

In the CT method we take the following separable expansion for the input kernel: 

5(2^:2/) = XI '^'"'(^<)^'(^>)' 2;|<| = {™ } (a;, y) (8) 



where 7;'s are expansion coefficients (not required to be determined), m;'s and wj's are defined by Eq. ([4]). 
The summation set S consists of the physical angular momenta {I}. 



^The requirements for the introduced input and output kernel quantities can be found in Ref[5]. We note that the 
Cox-Thompson expansions meet these criteria. 
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The requirements for the input kernel are i) to satisfy a differential equation (with a specific boundary 
condition) [5] and ii) to be uniquely solvable by the GL-type integral equation ([7]) at hand. While the first 
condition can be easily checked the second one is proved only for continuous transformation kernels [Sj. 
But this is not a restriction since K{x,y) is a twice continuously differentiable function. 

To solve the GL equation ([7|), Cox and Thompson[T] introduced a separable ansatz for the transfor- 
mation kernel, in the form of a sum over an artificial angular momentum space L G T 

K{x,y) = Y,AL{x)uL{y) (9) 

LeT 

where A^'s are unknown expansion functions and T is interpreted as a set of unknown shifted angular 
momenta to be determined under the constraint that it has the same number N of different elements as 
S does, i.e., iV = l^] = |T| and 5 n T = 0. 

By inserting the ansatze ^ and ^ into the GL equation (O and using the linear independence of the 
Riccati-Bessel functions, one obtains equations for the determination of the expansion functions Al{x) 
as follows 

Here, the only unknowns are the set T with elements L, and W denotes the Wronskian defined by 
W[a,h] = ab' -a'b. 

It can also be proved that Eq. (jlOp is uniquely solvable for continuous expansion functions, therefore 
the ansatz for the transformation kernel (O does not restrict the solution of the GL-type integral equation. 

If the set T is given one only has to solve a system of linear equations and perform a derivation in 
order to create the inverse potential Consequently the non-trivial task is to obtain T. 

In order to determine the set T one makes use of the Povzner-Levitan representation of the radial 
scattering wave function ([S]) whose asymptotic form, containing the input phase shift data, takes the 
form 

Bi sinix -l^ + 6i)= sm{x ' ) ' E ^^(^) i^i ""+1) I l^^lI 1) ' ^ ^ ^^^^ 

Here, Bi 's are normalization constants and we have defined the asymptotic expansion functions A'^ (x) = 
Al {x — >■ oo) which can be calculated from the asymptotic version of Eq. (|10[) which is given by 

Using the last two equations, (ITlT) and (TT^ . one can easily derive the following equations for the deter- 
mination of the unknown Us from the input phase shifts, 5i's[2j|3]: 

Si = \^ or tan(^0 = ,^^;t.+ ^V-V ^^^^ 
where Si = e^**' and the "shifted" reactance matrix elements are defined as 

/C±= [^^sin]^L[AC]L^e±'^('-'')"/^ leS, (14) 

with 

M,i„ 1 1 \ sin((/-L)f) ] leSL^T {V^\ 

M_ /,^-L(L + !)-?(/ + 1)1 cos((;-L)f) /' 'e^^eJ. (15) 

Clearly, the appearance of inverse of matrix Afcos makes solutions of Eq. (|13l) for the unknowns L 
especially hard. In practice this means that finding its solution is not an easy task for any nonlinear 
solvers. 

But, for instance with the deployment of sophisticated numerical methods, that can be applied to such 
ill-conditioned scenarios, one may succeed to solve either one or botlH of the highly nonlinear equations 
(|13p and thus find the set T. From that, as mentioned before, the CT inverse potential q(x) can easily 
be obtained by employing Eqs. (ITUl) . and dH])- 



^Despite the two sets of equations are equivalent with each other, for practical scenarios they are very different regarding 
their degree of ill-condition. 
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3 Semi-analytic solution 



In this section we present simplifications to Eqs. which can be used if only even (odd) partial 

waves are arising during the collision. The simplified equations can be employed also to construct a fair 
approximation to treat generic scattering problems. These approximations will be discussed in a later 
section. 

By solving Eq. (|12p explicitly without matrix inversion we will arrive to a simple set of nonlinear 
equations for the set T. First we will show that Eq. (fT2|) can be solved uniquely. This uniqueness is a 
consequence of the following statement. 

Lemma 3.1 For Xi,yj G C Xi ^ yj G / (\I\ = N) the system 

E^^ = 0' (16) 

has the unique solution Oi = \li ^ I . 

In the proof of the uniqueness we consider two solution sets and which satisfy Eq. 

(fT2)) . For the differences {A\{x) — we have a set of equations of the type (IT6l) . Hence because of 

Lemma 1 the difference between two sets of solutions is zero which means the uniqueness of the solution 
ofEq. (US]). 

Now let us differentiate Eq. twice with respect to the variable x. Then by employing the 

uniqueness result we arrive at the following equation 



dx^ 

which has a periodic solution as 



-Alix), (17) 



A^(a;) = flL cos(x) + &L sin(a;). (18) 

Now, by inserting this solution ([T51) into Eq. (|12l) and taking into account the independence of the sine 
and cosine functions, one gets the following two sets of equations 

yfaLl COS ((?-£)§) fcos(;f)1 



Consider two distinct cases of input data. In the first case let S consist of only even numbers 
(henceforth let us denote such a set with 5c) and in the second only of odd numbers (^o). In these two 
distinct cases, instead of (1191) . we have two kinds of equations (notice the simplification in comparison to 

my- 



(if) _ r 1 



^^^^ , bL \ L{L + l)-l{l + l) \ 

and 



L J L{L + 1) - l{l + 1) 



where \Tc\ = \Sc\, \To\ — \So\ and To n 5*0 = 0, To n = 0. In order to solve these sets of equations we 
need another statement similar to the one for the sets of equations of type (fTO)) . 

Lemma 3.2 For Xi,yj € C Xi ^ yj \/i,j G / (\I\ — N) the system 

E^^--l' (22) 



has the unique solution 



a, = n^g/("^ k e L (23) 
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Employing Lemma 1 and 2 for Eqs. ((20|) and (|2T]) one can arrive at the solutions 

UiesML + l)-l{l + l)) 1 
ilL'en\{L}iHL + 1) - + 1)) cos (if) 

and 

llL'eTo\{L}(^(-^ + 1) - ^ + 1)) sm (Lf ) 

respectively. 

Now, by using the explicit expressions ([M)) in Eqs. (IT51) and pT|) . one obtains finally the simplyfied 
'semi-analitic' equations to the CT method as 

valid for the case of even Z's. Similarly, using Eqs. ([25]) we get the semi-analytic equations to the CT 
method as 

which are valid in the case of odd Vs. These equations determine the unknown set To or To of shifted 
angular momenta 

Notice the simplified structure of the nonlinear equations and (j27p . compared to Eqs. (jl3D . 
While Eqs. (jl3p contain an explicit matrix inversion of a matrix involving the unknowns of shifted 
angular momenta, L's, formulae (1261) and (j27p do not require such a nonlinear operation (hence we 
use the term semi-analytic equations for it). They 'only' contain products and the tangent (cotangent) 
operations and are thus presumably easier to be solved for the sets Tq or To, if the respective input phase 
shifts are given. 

Finding the sets To or To, the corresponding potentials qdx) or qo{x) can be obtained similarly as in 
the general case, by employing Eqs. ([T0|. ([9]), and (jH). 



4 Applications 

As mentioned in the introduction the method discussed in section 3 can be applied for identical particle 
scattering of bosonic and fermionic type. It was checked that with the new equations such problems 
can be handled easily and without the deployment of sophisticated numerical methods, like simulated 
annealing [7]. 

Here we introduce another application of our method. Regarding the generic inverse scattering prob- 
lem consider the decomposition 

S — Sq lj 5*0 

where Sc and So contains, respectively, the even and odd elements of S. For the sets Sc and 6*0 one can 
perform the semi-analytic process. The resulting potentials are qc{x) and qo{x), respectively. 

Let us treat the general CT equations (|13l) as a set of coupled equations for the elements of T. Now 
if one disregards the coupling between the equations indexed by even and odd physical angular momenta 
then one gets an approximate method. It is obvious that in this approximation Tg U To is obtained for 
the set of the shifted angular momenta but let us use the sets Tc and To separately in the following. 

Considering, moreover, the simple additional prescription in equation ([9]), we conclude that it is 
reasonable to approximate the potential q{x) as 

ga(a;) = qdx) + qo{x). 

Since the solution of the semi-analytic equations is expected easier to be solved than the ones arising 
in the original procedure, the proposed approximation could be used for a larger class of input data. 
Below we demonstrate the applicability of the proposed approximation through inversion of experimental 
phase shift data. 

Chen and Thornow[^ have performed a comprehensive study of n scattering by ^^C target nucleus 
in the energy region 7 MeV < E}^^ < 24 MeV. They have derived 88 sets of complex-valued phase 

•^One can also determine the normalization constants. In both cases the result is _B; = [cos((5j)]~^. 
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shifts, one of which we use in our calculation here. Because of the spin-orbit coupling, each partial wave 
provides two phase shifts, and . In case of a weak spin-orbit coupling the combined phase shifts 
Si = [{I + l)S'^ + ld^]/{2l + 1) are characteristic of the underlying central potential^. 

Our input data will be the set of phase shifts at energy E^^^ = 12 MeV. In Table 1 we listed the 
quantities (5°"^ = KeSi, the real parts of the combined phase shifts, and also r/°"^ = \ exp{2idi)\, the 
elasticities. The results of the general CT procedure and the approximate method are visualized in Fig. 
1. 

Table 1: Original phase shifts 5°"^ and elasticities taken from the experimental phase shift data [see 
the text] for n scattering by ^^C nucleus at the energy E]^^ = 12 MeV. Both the real and imaginary parts 
of the shifted angular momenta L corresponding to the general (g) and the approximate (a) method, and 
the original phase shifts with the absolute differences between them and the calculated ones are listed. 
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A? = |<5?-5r|and - = |r;^-^r|. 



Figure 1: Inverse potentials V{r) obtained from input phase shifts (5°"^ and elasticities ri°"^ [see Table 
1] as a function of the radial distance r at energy E}^^ = 12 MeV (E','; "^- = 11.08 MeV, k = 0.766 fm"!). 
Potential components obtained by the generic method and the approximate one are denoted by solid and 
dashed curves, respectively. 
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To study the applicability of the approximate procedure one can compare the curves corresponding to 
the different procedures as shown in Fig. 1. We assume that the general method gives the correct result. 
As it can be seen the results obtained by employing the semi-analytic equations give a fair approximation. 
Indeed, if one calculates back the phase shifts and elasticities from the inverse potentials Vg and 14, we see 
that the reproduction is the best for the CT potential. Such observations can be made also by studying 
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Table 1 where we have hsted the differences A; and between the original phase shifts 5°"^ and ij""^ 
and the re-calculated ones provided by the inverse potentials of the two different methods. From this 
and other examples studied so far we conclude that the proposed approximation to the CT method can 
be used for a global orientation about the nature of the underlying interaction. 

5 Summary 

A simplified version of the main equation (ITS]) of the Cox- Thompson inverse scattering method at fixed 
energy was introduced. The simplification means that unlike the generic equation the new formulae do 
not involve the inversion of a matrix which contains the variables to be determined. Consequently the 
new equations (|26|) and (|27)) are easier to be solved by usual nonlinear solvers (such as that based, e.g., 
on the Newton- Raphson procedure [3]). However this simplification is only possible in cases when the 
partial waves arising in the scattering event are of the same parity. Thus the new equations are valid 
for the scattering of identical bosons or fermions. Because, in practice, there are also general collision 
experiments, we introduced an approximate procedure which is based on the fact that the potential is 
expressed as a sum for the shifted angular momenta. By an example taken from nuclear physics we 
demonstrated the applicability of the approximate method. 

We note that a general semi- analytical method which can be applied for all physical partial angular 
momenta is not yet known. Its development would be a promising field of future study. 
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